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THEORY  AND  GROUP  VELOCITY  OF 
ULTRASHORT,  TIGHTLY-FOCUSED  LASER  PULSES 

1.  INTRODUCTION 

Much  recent  interest  has  arisen  in  the  study  of  ultrashort,  high-intensity  laser  pulses. 
Numerous  experiments  are  underway  on  applications  of  ultra-intense  pulses,  including 
particle  acceleration  [1-3],  x-ray  sources  [4-7],  laser  fusion  [8]  and  the  study  of  ultrafast 
phenomena  [9].  This  activity  is  largely  a  result  of  the  development  of  compact  solid-state 
lasers  based  on  chirped  pulse  amphfication  [10],  which  are  capable  of  producing  ultrahigh 
intensities  (  >  10^®  W/cm^)  and  ultrashort  pulse  lengths  (  <  100  fs).  In  order  to  analyze 
short  pulse  phenomena  and  applications,  it  is  necessary  to  have  accurate  analytic  expres¬ 
sions  describing  the  properties  of  ultrashort,  tightly  focused  laser  pulses.  For  example,  it 
may  be  possible  to  focus  an  intense,  short  laser  pulse  onto  a  group  of  electrons  (or  a  diffuse 
plasma  or  gas)  such  that  the  electrons  are  picked  up  and  accelerated  by  the  laser  piilse. 
Electron  acceleration  can  be  a  result  of  the  axial  ponderomotive  force  associated  with  the 
fast  rise  in  the  laser  intensity.  Whether  or  not  the  electrons  are  trapped  and  efficiently 
accelerated  by  the  laser  pulse  depends  on  the  pulse  profile  and  is  a  sensitive  function  of 
the  pulse  group  velocity.  Accurate  expressions  for  the  pulse  envelope  evolution  and,  in 
particular,  the  group  velocity  are  necessary  to  evaluate  such  processes. 

A  common  approach  to  the  study  of  laser  pulses  is  to  make  use  of  the  paraxial  wave 
equation,  for  which  there  are  well-known  analytical  solutions  [11-13].  Strictly  speaking, 
these  solutions  are  valid  only  for  long  laser  beams:  L  Zr,  where  L  is  the  characteristic 
length  for  axial  variations  in  the  pulse  envelope  (typically  on  the  order  of  the  pulse  length), 
Zr  =  'Krl/X  is  the  Rayleigh  length  (or  diffraction  length),  ro  is  the  minimum  spotsize  of 
the  pulse  at  focus  and  A  is  the  laser  wavelength.  Frequently,  however,  the  paraxial  solutions 
are  modified  and  used  to  model  short  laser  pulses.  This  can  be  done  by  multiplying  the 
paraxial  solutions  by  an  axial  profile  function  of  the  form  f{z  —  Vet),  where  z  is  the  distance 
along  the  axis  of  propagation,  Ve  is  the  velocity  of  the  envelope  and  f/\df/dz\  ~  L.  These 
modified  paraxial  solutions  may  not  be  accurate  for  very  short  laser  pulses  {L  <  Zr). 

Furthermore,  it  is  not  clear  what  value  of  Ve  to  use  in  the  modified  paraxial  solutions. 
It  is  also  not  clear  how  the  envelope  velocity  Ve  is  related  to  the  group  velocity  of  the  pulse 
centroid  Vg.  Accurate  expressions  for  Vg  are  not  known.  Intuitively,  a  laser  pulse  passing 
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through  focus  has  Vg  <  c.  In  the  paraxial  approximation,  the  laser  spotsize  evolves 
according  to  r®  =  +  This  imphes  that  the  laser  pulse  photons  are  travehng 

at  the  diffraction  angle  dj,  =  tq/Zr  with  respect  to  the  z  axis.  The  axial  component  of 
the  photon  velocity  is  then  Wg  —  ccosdd  —  c(l  -  B‘d/2).  A  similar  result  is  obtained  by 
analyzing  the  phase  of  the  paraxial  solutions.  In  either  case,  \  —  The  results 

below  indicate  that  this  value  differs  from  the  group  velocity  of  the  pulse  centroid  Vg  by  a 
factor  of  two. 

In  this  paper,  analytic  solutions  to  the  wave  equation  [see  Eq.  (1)  below]  describ¬ 
ing  ultrashort,  tightly  focused  laser  pulses  are  found  to  arbitrary  order  in  the  parameter 
X/2irL  <  1.  In  addition,  exact  expressions  which  describe  several  properties  of  short  pulses 
are  derived,  such  as  global  conservation,  pulse  length  evolution,  centroid  motion  and  the 
pulse  group  velocity.  For  example,  to  leading  order,  the  laser  spotsize  evolves  according  to 
fs  =  +  where  rj  =  {z  +  ct)/2  and  the  initial  conditions  are  chosen  such  that 

the  minimum  spotsize  tq  occurs  at  the  focal  point  z  =  0  and  t  =  0.  Hence,  for  a  fixed  time, 
the  spotsize  varies  throughout  the  pulse.  Furthermore,  the  pulse  group  velocity  in  vacuum, 
based  on  the  motion  of  the  pulse  centroid,  is  given  by  1  —  Vg/c  ~  B'^/A,  to  leading  order. 
Detailed  comparisons  are  made  between  the  analytic  solutions  and  numerical  solutions  to 
the  full  wave  equation.  The  short  pulse  solutions  discussed  in  this  paper  are  hmited  to 
profiles  which  are  Gaussian  in  the  transverse  direction,  i.e.,  the  laser  field  amplitude  is  pro¬ 
portional  to  exp(— r^/r^),  where  r  is  the  radial  coordinate.  Generalization  of  the  results 
below  to  describe  higher  order  Hermite-Gaussian  transverse  profiles  is  straightforward. 

The  wave  equation  describing  the  three-dimensional  (3D)  evolution  of  of  a  laser  pulse 
in  a  fuUy  ionized  plasma  is  given  by 

where  a  =  eAj./meC^  is  the  normahzed  transverse  vector  potential  of  the  laser  pulse  (the 
transverse  electric  field  is  given  by  eEj_  =  —rriecda/dt)  and  kp  is  the  effective  plasma 
wavenumber.  In  general,  kp  is  a  highly  nonhnear  function  of  a  [3].  The  present  study, 
however,  will  be  Hmited  to  a  Hnear  plasma  response  where  kp  =  Wp/c  is  constant,  which 
is  vaUd  when  1,  where  Up  =  {4-7re^nofmey/^  is  the  plasma  firequency  and  no  is 
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the  ambient  electron  plasma  density  (assumed  to  be  uniform).  Here,  is  related  to 
the  intensity  of  a  hnearly  polarized  laser  field  by  ~  0.72  x  where  A  is  the 

laser  wavelength  in  /rm  and  I  is  the  laser  intensity  in  W/cm^.  In  the  <C  1  limit  (i.e., 
I  10^®  W/cm^  for  A  ~  1  nm.),  the  kp  term  models  the  dispersive  effects  of  the  plasma 
electrons.  The  vacuum  limit  corresponds  to  kp  =  0.  Equation  (1)  indicates  that  the  initial 
polarization  of  the  laser  field  remains  unchanged  as  it  propagates.  The  axial  component 
of  the  field  can  be  found  by  requiring  V  •  E  =  0. 

The  majority  of  previous  theoretical  studies  of  tightly  focused  pulse  propagation  has 
been  concerned  with  solutions  to  the  paraxial  wave  equation  [11-13].  Paraxial  solutions 
are  discussed  in  detail  in  Sec.  2.  I.P.  Christov  [14]  found  analytical  solutions  describing 
pulse  evolution  which  are  vahd  only  far  from  focus,  z  ^  Zr.  Furthermore,  the  group 
velocity  of  the  pulse  was  not  considered.  More  recently,  Horvath  and  Bor  [15]  described 
pulse  distortions  which  arise  when  a  pulse  is  focused  by  a  lens  with  a  frequency  dependent 
index  of  refraction.  These  distortions  are  due  to  the  longitudinal  chromatic  aberration 
of  the  lens,  i.e.,  the  focal  length  of  the  lens  is  a  function  of  frequency.  A  short  pulse  of 
length  L  has  a  finite  bandwidth  Auj  ~  1/X.  Hence,  the  different  firequency  components 
focus  along  different  paths  which  distorts  the  pulse  profile.  Pulse  distortions  due  to  lens 
aberrations  will  be  neglected  in  the  following. 

The  remainder  of  this  paper  is  organized  as  follows.  Section  2  discusses  solutions  to 
the  paraxial  wave  equation  and  their  shortcomings.  Section  3  describes  the  formalism  for 
obtaining  short  pulse  solutions  to  the  full  wave  equation.  Use  is  made  of  the  independent 
variables  ^  =  z  —  ct  and  rj  =  {z  +  ct)/2.  The  wave  equation  for  the  pulse  envelope  a(r,  rf,  C) 
is  solved  by  taking  a  Fourier  transform  with  respect  to  C)  where  a  is  given  by  the  real  part 
of  aexp(zfcoC)6x5  with  fco  =  27r/A.  Several  exact  properties  of  short  pulses  are  derived, 
including  global  conservation,  pulse  centroid  motion  and  group  velocity,  and  pulse  length 
evolution.  Section  4  presents  approximate  short  pulse  solutions  based  on  the  expansion 
parameter  1/koL  <  1.  The  zeroth,  first,  second  and  fourth  order  solutions  are  highhghted. 
Comparisons  between  the  analytic  solutions  and  numerical  solutions  to  the  wave  equation 
are  presented  in  Sec.  5.  The  paper  concludes  with  a  discussion  in  Sec.  6. 
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2.  PARAXIAL  SOLUTIONS 


The  mostly  commonly  used  expressions  describing  the  evolution  of  laser  pulses  are 
the  solutions  to  the  paraxial  wave  equation  [11-13].  The  paraxial  wave  equation  can  be 
obtained  as  follows.  Consider  a  long  laser  pulse,  i.e.,  a  laser  “beam”,  of  the  form  a  == 
a{r,  z)  exp  [ikQ{z  -  ct)]  e^,  where  d(r,  z)  is  the  laser  envelope  and  cko  is  the  laser  frequency. 
The  wave  equation  describing  the  laser  envelope  is  given  by 

^]_  +  2iko—  +  -^-kp^a  =  0.  (2a) 

The  paraxial  approximation  involves  neglecting  the  d'^/dz^  term,  which  assumes  a  suffi¬ 
ciently  slowly  varying  envelope,  \dal dz\  <^\kQa\.  The  paraxial  wave  equation  is 

(vi  +  2iA:o^-A:2^a  =  0.  (26) 

Solutions  to  the  paraxial  wave  equation  are  well-known  in  terms  of  Hermite-Gaussian 
modes  [11-13].  For  example,  the  lowest  order  Gaussian  mode  is  given  by  d  =  ap{r,kQ,z), 
where  dp  =  uq  exp{ipp)  and 

i^p{r,ko,z)  =  -^ln(l  +  Q!2)  -  a ^  -  jk^rla^.  (3) 

In  Eq.  (3),  =  z/Zr,  Zr  =  fcoro/2  is  the  Rayleigh  length,  gq  is  the  peak  amphtude 

at  focus,  and  tq  is  the  minimum  laser  spotsize  at  focus.  For  simplicity,  the  focal  point  is 
chosen  to  be  at  z  =  0.  The  properties  of  the  paraxial  solutions  are  well-known  [11-13].  For 
example,  |dp|  =  (aoro/rj  exp(-r2/r2),  where  r*  =  ro(l  -f-  z^jZ]^'^!'^  is  the  laser  spotsize. 
The  paraxial  solution  Op  =  dp  exp  [ikQ{z  —  ct)]  can  be  written  as 

ap{r,z,t)  =  {aQro/rs)exp{-r^ /r^  +  i(t)p),  (4a) 

<j)p{r,  z,  t)  =  ko{z  -  ct)  -I-  OL^r'^Irl  -  tan“^  -  k^r^ajA,  (46) 


where  <pp  —  [ipp]i  +  ko{z  -  ct)  is  the  total  phase  (the  subscript  i  denotes  the  imaginary 
part).  The  effective  pulse  frequency  u;  and  axial  wavenumber  kz  can  be  defined  in  terms 
of  the  total  phase,  uj  =  -d(j)p/dt  and  kz  =  d^p/dz.  This  gives  uj  =  cko  and 
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The  paraxial  approximation  implicitly  assumes  k^/kQ  1  and  kQvl  ^  1. 

It  is  reasonable  to  assume  that  if  the  laser  pulse  length  L  is  sufficiently  long,  L  ^ 
Zr  ^  A,  where  A  =  27r/A:o  is  the  laser  wavelength,  then  the  pulse  can  be  adequately 
described  by  a  paraxial  solution  of  the  form 


a(r,  2,  t)  =  f{z  -  Vet)ap{r,  z,  t) 


(6) 


Here,  the  function  /  describes  the  axial  profile  of  the  pulse  envelope,  which  is  assumed  to 
travel  at  the  envelope  velocity  Vg,  and  f/\df/dz\  ~  L.  It  is  not  clear,  however,  what  value 
of  Vg  to  use  in  the  axial  profile  f{z  —  Vgt).  Within  the  paraxial  approximation,  the  “group 
velocity”  of  the  pulse  envelope  Vg  can  be  estimated  from  the  total  phase  (f)p  of  the  laser 
field.  In  terms  of  the  frequency  u  =  cko  and  axial  wavenumber  Eq.  (5),  the  effective 
phase  velocity  Vp  and  “group”  velocities  Vg  are  given  by  Vp  =  w/kz  and  Vg  —  {dkzldu})~^ . 
Specifically,  Vg/c  =  (1  +  ziil  —  €p,  or 


Ve{r,  z) 
c 


(1  —  6q:^  +  af) 
^  (l  +  a2)3  _ 


(7) 


In  particular,  in  vacuum  kp  =  Q  and  at  r  —  z  =  0,  Vg/c  ~  1  —  2/A:oro.  Note  however, 
that  this  expression  for  Vg  is  problematic,  since  it  imphes  that  there  are  regions  in  {r,z) 
for  which  Cp  <  0  and  Vg  >  c.  To  provide  an  accurate  description  of  the  group  velocity,  it  is 
necessary  to  consider  the  motion  of  the  pulse  centroid,  as  is  discussed  below. 


3.  EXACT  DESCRIPTION  OF  SHORT  PULSES 

To  describe  the  behavior  of  short  pulses,  it  is  convenient  to  introduce  the  variables 
^  =  z  —  ct  and  rj  =  {z  +  ct) /2.  In  these  variables,  the  wave  equation  becomes 

Short  pulse  solutions  will  be  sought  of  the  form  a  =  [dexp(iA:oC)  +  c-c.]  ea./2,  where 
C»  v)  is  the  laser  envelope,  fco  is  a  constant  and  represents  the  fundamental  laser  pulse 
wavenumber,  is  a  unit  vector  in  the  direction  of  the  polarization  and  c.c.  denotes  the 
complex  conjugate.  Using  this  form,  the  laser  pulse  propagates  in  the  positive  z  direction. 
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C  is  constant  for  a  point  moving  at  the  speed  of  light  and  is  a  measure  of  the  relative  axial 
distance  within  the  laser  pulse,  and  77  is  the  time-like  coordinate  which  represents  the  axial 
distance  traveled  for  a  point  moving  at  the  speed  of  hght.  The  wave  equation  describing 
the  evolution  of  the  envelope  d  is  given  by 

Taking  the  Fourier  transform  in  C  gives 

Vi  +  2i(i„  +  k)^-  kl 

where 

1 

is  the  Fourier  transform  of  a. 

Notice  that  Eq.  (9b)  is  identical  to  the  paraxial  wave  equation,  Eq.  (2b),  with  z  rj 
and  A:o  ^  A:o  +  k.  Hence,  explicit  solutions  to  Eq.  (9b)  can  be  found  based  on  the  Hermite- 
Gaussian  solutions  of  the  paraxial  wave  equation  [11-13].  For  example,  assuming  the  lowest 
order  Gaussian  mode,  the  Fourier  transform  of  the  laser  envelope  is  given  by 


dfc(r,fc,77) 


/: 


dC  exp{-ikC)d{r,  C,  v) 


(9c) 


dfc  —  0, 


(96) 


a{r,C,v)  =  0. 


(9a) 


ak{r,  k,  Tj)  =  fkdpk{r,  ko  +  k,rj)  =  aofk  expipkir,  ko  -f  fc,  77),  (10a) 

2/2 

^kir,ko  +  k,v)  =  -^ln(l  +  a^)-  -itan'^a^  -  ^k^r^ak,  (106) 

2  1  +  lOLk  4 

ak  =  vl^Rk,  (10c) 

ZRk  =  (ko  +  k)rl/2.  (lOd) 


Note  that  dfc(77  =  0)  =  /^ao  exp(-r^/ro).  Hence,  fk  is  the  Fourier  transform  of  the  initial 
(77  =  0)  axial  envelope  profile  /(C)-  For  convenience,  the  focal  point  is  chosen  to  occur  at 
77  =  0. 

Equation  (10)  is  a  vahd  solution  for  the  lowest  order  Gaussian  mode  to  Eq.  (9b)  for  all 
k  except  k  =  -ko-  Notice  that  k  =  -ko  corresponds  to  a  Fourier  mode  component  which 
is  axially  uniform,  i.e.,  da/dC,  0,  which  is  physically  uninteresting.  When  k  =  -ko,  the 
d/dr]  operator  in  Eq.  (9b)  and  the  quantity  ock  are  singular  and  the  solution  Eq.  (10)  breaks 
down.  These  singularities  can  be  avoided  by  choosing  a  distribution  fk  which  does  not 
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contain  components  at  k  =  —ko.  For  analytical  convenience  in  the  following  calculations, 
an  axial  k  spectrum  will  be  used  of  the  form  fk  =  (1  +  k/ko)fGk,  which  corresponds  to 
/  =  (1  —  ikQ^d/dC)fG,  where  Jg  =  exp(— represents  a  Gaussian  axial  pulse  profile 
with  a  pulse  length  L.  Specifically,  the  following  forms  for  f{()  and  fk  will  be  assumed. 


expl-;p 


fk 


k 


e 
2 


(11a) 

(116) 

(12) 

where  Iq  is  the  intensity  of  the  pulse  center  at  focus  (at  focus,  the  pulse  center  is  given 
by  C  =  z  —  ct  =  0  and  r  =  0).  This  represents  a  slightly  distorted  Gaussian  profile,  since 
1/k^L^  <  1. 


Equation  (11a)  implies  an  intensity  profile  I  ~  |ap  at  focus  (rj  =  0)  given  by 

I  =  Io 


r, ,  4(z-ct)2] 

1 

to 

1 

a 

to 

to 

to 

kU‘  \ 

exp 

1 

CNO 

1  ^ 

1 

cs 

_ 1 

3.1.  Global  Conservation 


Based  on  the  wave  equation  for  the  pulse  envelope  a{r,(,v)y  Eq.  (9a),  it  is  straight¬ 
forward  to  show  that  the  quantity 


W  = 


(13) 


is  an  exact  constant  of  the  motion,  i.e.,  dW/dr]  =  0.  Physically,  W  may  be  interpreted 
as  the  leading  order  contribution  to  the  total  pulse  energy  at  a  given  r],  i.e.,  the  energy 
density  integrated  over  the  transverse  and  axial  coordinates.  In  normahzed  units,  the 
leading  order  components  of  the  the  transverse  electric  and  magnetic  fields  are  Ex  = 
c~^da/dt  —  (5/5C  “  2^/9t])o-  and  By  ~  da/dz  =  {d/dC,  +  |5/5ry)a.  Hence, 

-f  ~  |5a/5Cp  +  |5a/977p/4  ~  |5o/5CP 

which  assumes  ^  \daldx\^. 

Letting  h  =  (ifco  +  d/dOa  and  letting  hk  =  i{ko  +  k)ak  denote  the  Fourier  transform 
of  6,  W  can  be  evaluated  using  the  identity  f  dCl^P  =  J  dk\bk\^.  One  finds 

W  =  ^  r  dkiko  +  k)^\fk\\  (14) 
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where  use  has  been  made  of  the  fact  that  /  drr\apk\^  =  aoro/4.  For  the  specific  form  of 
fk  given  by  Eq.  (11), 

W  =  ^Ulrlkl  (l  +  ^  +  ^)  ,  (16) 

where,  typically,  —  (27rL/A)^  S>  1 

3.2.  Pulse  Centroid  and  Group  Velocity 

Since  |6p  is  a  globally  conserved  quantity,  it  can  be  used  as  a  weight  function  for 
defining  the  centroid  C  of  the  laser  pulse,  i.e., 

CM  =  (0  =  r  drr  r  dCC\h\VW.  (16) 

JQ  J  —  oo 

The  centroid  can  be  calculated  explicitly  using  the  relation  f  dCgb  =  f  dkgkbl,  where 
g  =  (b  with  its  transform  gk  =  idbkfdk  and  b*  represents  the  complex  conjugate  of  b. 
Assuming  /fe  to  be  real  gives 

1  f°°  r°°  f) 

(  =  -—J  drr  J  dk{ko  +  ky\dk\'^—lia{ipk) 

For  the  specific  form  of  fk  given  by  Eq.  (11), 

_  v{l  +  kyo/2){l  +  1/klL^)  ^  _ 

^  jfc2r2(l  +  6/fc2L2  +  3/fc^L4)  ^ 

Note  that  if  a  pure  Gaussian  axial  profile  fc  =  exp(-C^/I'^)  was  used  in  Eqs.  (14) 
and  (17)  instead  of  the  modified  Gaussian  profile  of  Eq.  (11a),  this  would  result  in  corrects 
of  order  l/k^L^  or  higher.  Specifically,  the  factor  (1  +  G/k^L^  +  S/k^L^^)  appearing  in  Eq. 
(15)  and  would  be  replaced  by  (1  +  2/k^L^)  and  Eq.  (18a)  would  become 


(18a) 


77(1  +  fcgrg/2) 


(185) 


klr^{l  +  2/k^oL^y 

The  evolution  of  the  pulse  centroid  C  can  provide  a  definition  for  the  group  velocity 
of  the  pulse,  where  Vg  =  dz/dt  with  ^  =  z  —  ct.  Setting  z  =  z  in  Eq.  (18a)  and  letting 


r]  —  (z  +  ct)/2  gives 


Vg/c  =  (1  -  e/2)/(l  +  e/2). 
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where  e  =  —dC/dr]  can  be  determined  trivially  from  Eq.  (18a).  In  the  hmits  e  <C  1  and 

»  1, 

1  /  k'^rl  \ 

(20) 


Vg/c  ~  1  -  e  ~  1  -  I  1  +  — 


0 


klL-^ 


The  relativistic  factor  7g  associated  with  Vg  vs.  'fg  =  (1  —  u^/c^)  ~  1/V^,  assuming 

e  <  1. 

3.3.  Pulse  Length 

In  a  similar  fashion,  the  evolution  of  the  pulse  length  Lp  can  be  determined  by  defining 
^P  =  ((C^)  -  where 


poo  poo 

(C")  =  /  drr  da^iWW. 

»/  0  J  “oo 

Noting  that  f  d^g‘g  =  J  dkg^gk  and  assuming  fk  real  gives 


(21) 


2  2  roo 

_  “o"o  / 

4ir  J_ 


dk 


{ko  +  k) 


+  fk 

f2„2 


dk 


d'ipk\ 


2/5^^^  (,,ky.ky. 


+ 


/  {ko  +  k)'^rQ 

For  the  specific  form  of  fk  given  by  Eq.  (11), 


1  + 


+ 


(22) 


(C^)  =  ^ 
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k^L^  k^L^  fcgL2 


^  2  8 


6  3  \ 


-1 


In  the  hmit  k^L^  ^  1  and  fcprg  -C  1, 


r2 

Lic:i  — 

P  4 


1  + 


^2X2 


+ 


nVzl  (,  i 

k^L^  V  k^L^ 


(23) 


(24) 


Hence,  the  pulse  length  slowly  increases  as  the  pulse  moves  away  from  the  focal  point 
T]  =  0.  In  particular,  Lp{ri)/Lp{0)  ~  1  +  [8Z^A:oLp(0)] ,  assuming  k^L"^  >  1. 

4.  APPROXIMATE  SOLUTIONS 

The  laser  pulse  envelope  a(r,  (,  rj)  is  given  by  the  inverse  transform  of  Eq.  (10), 

1  f°° 

a{r,C,'n)  =  -i=  dkexTp{ikC)fkaoexp{ipk)-  (25) 

y  27r  J — oo 


9 


For  a  particular  choice  of  fk,  it  is  unlikely  that  exact  analytical  solutions  can  be  obtained 
for  the  inverse  transform  Eq.  (25).  Physically,  k  represents  the  transform  of  the  C  variations 
which  occur  in  the  laser  envelope.  Typically,  k  ~  1/L  and  k/ko  ~  l/koL  <C  1.  Hence,  the 
function  'ipkiko  +  k)  can  be  expanded  about  ko  and  approximate  expressions  for  the  inverse 
transform  can  be  obtained,  i.e.. 


where 


9V'  ,  1  ,_2 


V’fc(r,  ko  +  k,‘n)  =  ^(ko)  +  —k  +  - 

OKq  Z  t//VQ 


r  + 


V'(fco) 

dij) 

dko 

dkl 


-iln(l  +  a2) 


la 

ko 


— itan  ^a—^kpr^a, 
\  +  ia  4  P  ’ 


+ 


la 

Kq 


1  +  ZQ!  (l  +  za:)2  '  4 

2  +  ia  2r’^/ro 

(1  +  (l  +  io:)^ 


+ 


PP 


(26o) 


(266) 

(26c) 

(26d) 


and  a  =  'hIZr. 


4.1.  Zeroth  Order  Solution 


To  lowest  order,  tl)k  —  V'(^o),  independent  of  k.  The  inverse  transform  Eq.  (25)  gives 
the  zeroth  order  solution 


aS^\r,  C,  v)  =  aofiO  exp(V').  (27) 

The  zeroth  order  solution  is  similar  to  the  paraxial  solution,  Eq.  (6),  with  the  independent 
variable  z  replaced  by  77  =  (z  +  ct)/2.  Furthermore,  to  zeroth  order,  the  axial  profile  is 
=  /(z  —  ct),  i.e.,  the  envelope  propagates  with  Vg  =  c.  Notice  that  the  initial  {t  —  0) 
envelope  is  given  by  a{t  =  0)  =  00/(2;)  exp  [V’(z/2)],  in  contrast  to  the  paraxial  solution, 
for  which  dp(t  =  0)  =  aof{z)  exp  [i}{z)\. 

The  full  zeroth  order  solution  o(°)  =  exp(iA:oC)  can  be  written  as 


=:  koC  +  Oir'^/r^  —  tan  ^  a  —  kpr^a/A, 
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(28o) 

(286) 


where  +  fcoC  is  the  total  phase  of  the  zeroth  order  solution  and  =  ro(l  + 

is  the  spotsize.  The  effective  pulse  frequency  and  axial  wavenumber  of  the  zeroth 
order  solution  can  be  defined  by  0;^°^  =  —dcp^^^/dt  and  =  d(j)^^^/dz.  This  gives 


c.(°)  =  cfco<!l+-^  + 


4^0  L 


l-3(l-a2) 


=  koll- 


L 


(29a) 

(296) 


Notice  that  Eqs.  (29a)  and  (29b)  satisfy  the  local  dispersion  relation 


P  '  rp2 


(30) 


This  is  in  contrast  to  the  paraxial  solution,  in  which  the  above  dispersion  relation  is  only 
approximately  satisfied  for  lj  =  cko  and  k^  given  by  Eq.  (5). 


4.2.  First  Order  Solution 


To  first  order,  -ipk  —  ip{ko)  +  ^'fc,  where  the  prime  denotes  d/dko.  The  first  order 
solution  is 

_  _J_ 

=  Oo/(C  -  i^')  exp(V').  (31) 


C,  r}) 


/: 


dk  exp  [ifc(C  -  i'ij}')]  fkdo  exp(^) 


The  first  order  solution  includes  corrections  to  the  evolution  of  the  axial  profile,  i.e.. 
In  particular,  first  order  group  velocity  effects  are  included.  Intuitively,  the  local 
envelope  velocity  can  be  estimated  by  setting  d(C  -  iip')r/dt  =  0,  where  the  subscript  r  de¬ 
notes  the  real  part,  and  by  identifying  Ve  =  dz/dt.  This  gives  Ug/c  =  (1  -  ei/2)/(l  +  ei/2). 


i.e.. 


Ve{r,  v) 


c 

T-l: 


~  1  —  Cl  =  1 


l.2j.2 
”'0' 0 


(1  -  a^)  r\l  -  6a"  +  a^)  ,  k^r^ 


+ 


(32) 


(l-fa^)^  r§(l -I- 0:2)3  4 

where  €1  =  Z^^d'ipyda  (the  subscript  i  denotes  the  imaginary  part).  Alternatively,  a  local 
value  of  the  pulse  centroid  as  a  function  of  radius  can  be  defined  by  averaging  over 
the  axial  profile,  i.e.. 
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where  a  Gaussian  axial  profile  was  assumed,  /  =  exp(-C^/i^).  (If  the  quantity  \b\^  is  used 
as  a  weight  function  with  /  given  by  Eq.  (11),  an  identical  result  is  obtained  plus  corrections 
of  order  l/kgL^  or  higher.)  Equation  (33)  imphes  that  Ve/c  =  (1  —  ei/2)/(l  +  ei/2)  as 
before,  i.e.,  the  local  envelope  velocity  is  given  by  Eq.  (32). 

The  above  discussion  sheds  some  light  on  the  interpretation  of  the  paraxial  result  for 
the  envelope  velocity  given  by  Eq.  (7).  Notice  that  the  paraxial  result  is  identical  to  Eq. 
(32)  with  ri  -*  z  in  the  definition  of  a.  Both  Eqs.  (7)  and  (32)  exhibit  regions  in  (r,  a) 
for  which  Ve  >  c.  Hence,  neither  of  these  expressions  should  be  interpreted  as  the  pulse 
group  velocity.  Instead,  Eqs.  (7)  and  (32)  should  be  interpreted  as  the  velocity  of  the  local 
pulse  centroid  (C)c  defined  at  a  specific  radius.  Note  that  a  change  in  the  local  value  of 
(C)c  at  fixed  r  does  not  necessarily  correspond  to  an  axial  transport  of  energy.  Changes  in 
(C)c  can  reflect  distortions  in  the  pulse  profile  as  it  evolves.  Consider,  for  example,  a  pulse 
profile  which  is  initially  symmetric  and  centered  about  the  focal  point  such  that  (C)c  =  0 
at  all  T.  As  the  pulse  propagates,  the  front  of  the  pulse  expands  radially  (diffracts)  as  it 
moves  away  from  the  focal  point.  The  back  of  the  pulse,  which  is  initially  located  behind 
the  focal  point,  wiU  contract  radially  as  it  moves  towards  the  focal  point.  The  pulse  is 
now  asymmetric,  with  the  front  of  the  pulse  wider  than  the  back.  For  sufl&ciently  large 
r,  this  can  cause  the  local  centroid  (C)c  to  move  towards  the  front  of  the  pulse  resulting 
in  a  “local”  Ve{r,r})  >  c,  as  indicated  by  Eq.  (32).  However,  this  local  centroid  motion  is 
associated  with  the  distortions  of  the  pulse  profile  and  a  radial  transport  of  energy,  not 


with  an  axial  transport  of  energy. 

It  is  more  accurate  to  identify  the  group  velocity  of  the  pulse  with  the  motion  of  the 
actual  centroid,  C>  of  the  entire  pulse.  The  pulse  centroid  C  =  (C)  is  given  by  averaging 
over  both  the  radial  and  axial  profiles: 


r  _  f  drr  JdCC|/Pexp(2V>^)  _ _ ^ 


C  = 


^0^0 


1  + 


e 


P'O  \  _ 


=  -e2»7. 


(34) 


J  drr  f  dCl/Pexp(2V'r-) 

To  first  order,  Eq.  (34)  is  equivalent  to  Eq.  (16),  producing  the  leading  order  contribution 
to  the  exact  expression  given  by  Eq.  (18a).  The  group  velocity  of  the  centroid  is  then 


Vg/C 


(1-^2/2) 

(1  +  ^2/2) 


~  1  — 


1  + 


(35) 
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which  is  a  constant  independent  of  77. 

In  addition  to  group  velocity  effects,  the  first  order  solution  has  additional  terms 
which  affect  the  phase  of  the  laser  pulse.  For  example,  the  leading  order  correction  to 
the  total  phase  is  given  by  +  2C^'/L^,  where  0^°^  is  given  by  Eq.  (28b). 

This  additional  term  gives  contributions  to  the  effective  frequency  and  wavenumber  of 
the  pulse.  For  example,  the  effective  frequency,  defined  as  =  —d<f)^^ydt,  is  given  by 
where  0;^°^  is  given  by  Eq.  (29a)  and 

_  2  f 

cko  klL^  I  (1  + 

Notice  that  along  the  axis,  r  =  0,  the  frequency  varies  slightly  throughout  the  pulse, 
=  Su;(^^(C),  and  evolves  as  the  pulse  propagates  away  from  focus, 

4.3.  Second  Order  Solution 


2r^ 


Ca 


Znil  +  a^y  L 


1-  (l_a») 


(36) 


For  a  pulse  with  a  Gaussian  axial  profile,  analytical  expressions  can  be  obtained  for  a 
to  second  order.  Approximating  V’fc  =  V’  +  +  il)"k^/2  and  using  the  expression  for  fk 

given  by  Eq.  (11)  gives 


C,/?)  = 


ao 


1  + 


2i(C  -  #') 


^1  -  2r/L^  L  koL\l-2riL-^)\ 


exp 


1}) 


(l-2^'V-^2) 


•  (37) 


This  analytic  form  for  the  second  order  solution  will  be  compared  with  numerical  solutions 
to  the  wave  equation  in  Sec.  5. 


4.4.  Solution  to  Arbitrzury  Order 


Keeping  the  full  expansion  for  ■0*.  and  expanding  the  second  order  and  higher  terms 
in  the  exponential  gives  a  solution  to  arbitrary  order, 

aW(r,  C,v)  =  ^j°°  dk(l  +  S+^S^  +  ---  +  exp  [7/;  +  i{C  -  i^')k] ,  (38) 


where 


5  = 


1  ,.2 


2  dk^ 


r  + 


1  ^ 

^  n!  dk^  ■ 


(39) 


Equation  (38)  can  be  readily  evaluated.  For  example,  keeping  all  terms  up  to  order 


d^^\r,C,v)  =  0,0^^ 


r,  r 

a,"'  a"  ^ 

irf ' 

an 

^  2  dcp- 

+  6  ac^'^ 

[  24 

8  . 

- 1 

/(C  -  iV'')-  (40) 
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5.  COMPARISON  TO  NUMERICAL  SOLUTIONS 


The  wave  equation,  Eq.  (9a),  can  be  solved  numerically,  using  standard  techniques, 
to  obtain  a(r,  rj).  For  example,  consider  a  laser  pulse  at  the  point  of  focus  {rj  —  0)  with 

d(r,C,0)  =  a„(l  +  ^)  exp  [-  (^  +  ij)]  .  (41) 

The  evolution  of  this  laser  pulse  in  vacuum,  as  a  function  of  r],  can  be  estimated  using 
the  first  order  solution  Eq.  (31)  [with  /(^)  =  (1  +  2z^/A;oL^)  exp(— C^/X^)],  the  second 
order  solution  Eq.  (37),  or  by  numerically  solving  Eq.  (9a)  using  Eq.  (41)  as  the  initial 
condition.  The  result  for  the  second  order  solution  is  shown  in  Fig.  1,  where  the  real  part 
of  o  =  dexp(ifcoC)  is  plotted  versus  r  and  C  at  fixed  time  {t  =  rj  -  (/2)  for  (a)  t  =  0  and 
(b)  t  =  2Zii/c.  In  this  example,  A  =  2.5  //m  and  L  =  ro  =  b  /xm,  such  that  Zr  ~  31.4 
and  L/X  =  tq/X  =  2.  The  curvature  of  the  wave  fronts  as  the  laser  pulse  moves  away  from 
focus  is  clearly  observable  in  Fig.  1(b). 

The  numerical  solution  for  the  parameters  of  Fig.  1  is  compared  to  the  first  and  second 
order  solutions  in  Fig.  2.  Here,  the  percentage  difference  100(a(”)  —  a^”“’”))/Max(a("^)  is 
plotted  versus  r  and  C  at  77  =  Zr  for  (a)  the  first  order  solution  and  (b)  the  second 
order  solution  where  denotes  the  numerical  solution  and  Max(a(”))  denotes 

the  maximum  value  of  the  n-th  order  solution  The  numerical  solution  agrees  with  the 
approximate  solutions  to  better  than  1%.  Numerical  results  show  that  this  error  decreases 
as  L/X  increases. 

The  group  velocity  of  the  numerical  solution  can  be  determined  by  numerically  eval¬ 
uating  the  pulse  centroid  C,  Eq.  (16),  using  the  numerical  solution  in  the  integrand.  From 
Eq.  (19),  the  group  velocity  is  given  by  Vg/c  =  (1  —  e/2)/(l  -f  e/2),  where  the  quantity 
e  =  —d^/dr]  is  evaluated  numerically.  The  result  is  compared  to  the  analytical  solution, 
Eqs.  (18a)  and  (19),  in  Fig.  3,  where  1  —  Vg/c  is  plotted  versus  rj  for  each  case  using  the 
same  parameters  as  in  Fig.  1.  The  local  envelope  velocity  Ve  a.t  r  =  0  can  be  similarly 
determined  from  the  numerical  solution  by  integrating  only  over  C  on  the  right-hand  side 
of  Eq.  (16).  The  result  is  compared  to  Eq.  (32)  in  Fig.  3.  The  numerical  results  show 
excellent  agreement  with  theory. 
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6.  SUMMARY  AND  DISCUSSION 


New  solutions  to  the  wave  equation,  Eq.  (1),  have  been  derived  which  describe  an 
ultrashort,  tightly  focused  laser  pulse  propagating  through  a  plasma  with  a  linear  re¬ 
sponse  or  through  vacuum.  The  wave  equation  for  the  pulse  envelope  a{r,(^,rj),  where 
a  =  aexp(tA:oC)  +  c.c,  was  solved  making  use  of  the  independent  variables  C  =  z  —  ct  and 
77  =  {z  +  ct)/2  and  by  taking  a  Fourier  transform  with  respect  to  C-  Solutions  for  the  trans¬ 
form  of  the  envelope,  dfc(r,  ko  -I-  k,  77),  could  then  be  found  without  approximation  in  terms 
of  Hermite-Gaussian  modes.  The  solutions  are  analogous  to  solutions  of  the  paraxial  wave 
equation.  For  simphcity,  this  study  was  hmited  to  envelopes  in  which  the  radial  profile  is 
described  by  the  fundamental  Gaussian  mode,  i.e.,  |a|  ~  exp(— r^/7*^).  Generahzation  to 
higher  order  Hermite-Gaussian  modes  should  be  straightforward.  Letting  the  solution  to 
the  paraxial  wave  equation,  Eq.  (2b),  be  denoted  by  d  =  dp(r,  ko,z),  the  exact  solution  to 
the  full  wave  equation  for  d^,  Eq.  (9b),  is  given  by  d*,  =  /fcdp(r,  ko  -1-  A:,  77),  as  is  indicated 
by  Eq.  (10),  where  /*,  is  the  fourier  transform  of  the  axial  pulse  profile  at  77  =  0. 

Based  on  the  exact  solution  for  the  transform  of  the  envelope  d^,  several  exact  prop¬ 
erties  of  ultrashort,  tightly  focused  laser  pulses  were  determined.  General  expressions  were 
determined  for  arbitrary  axial  profiles,  fk,  and  specific  expressions  were  determined  for 
the  case  of  a  quasi-Gaussian  axial  profile  given  by  Eq.  (11). 

For  a  pulse  with  a  quasi-Gaussian  axial  profile,  the  quantity  W  =  J  drr  f  dC\b\^,  where 
b  =  {iko-\-d/dC)d,  was  identified  as  a  globally  conserved  quantity  of  the  full  wave  equation 
for  the  envelope,  Eq.  (9a),  i.e.,  dW/dr]  =  0.  The  quantity  W  is  roughly  equivalent  to  the 
total  pulse  energy  at  fixed  77.  Using  the  quantity  |&p  as  a  weight  function,  an  expression 
for  the  pulse  centroid  (  was  determined,  Eq.  (18a),  by  integrating  (  over  the  axial  and 
radial  coordinates.  The  motion  of  the  centroid  was  used  to  define  the  group  velocity  of 
the  pulse,  Vg  =  dz/dt,  where  (  =  z  -  ct.  The  resulting  value  for  Vg  is  given  by  Eq.  (19). 
By  averaging  over  the  axial  and  radial  coordinates  with  as  the  weight  function,  an 
expression  for  the  laser  pulse  length  Lp  was  determined,  Eq.  (24).  It  was  found  that  the 
pulse  length  slowly  increased  as  the  pulse  propagates. 

Approximate  solutions  to  the  laser  pulse  envelope  d(r,  C,  rj)  were  found  by  expanding 
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ak{r,ko  +  k,r})  about  ko  for  |/i:/fco|  1/koL  1  and  evaluating  the  inverse  transform. 
Approximate  solutions  for  a{r,C,T])  were  found  to  arbitrary  order  in  \k/ko\.  The  zeroth 
order  solution,  Eq.  (27),  gave  =  /(C)flp(’')  ^o,  ^),  where  ap{r,ko,z)  is  the  solution  to 
the  paraxial  wave  equation.  This  indicates  that  in  order  to  describe  a  short  pulse  with  the 
paraxial  solution,  the  variable  2  should  be  replaced  by  the  correct  propagation  variable 
77  =  (z  +  ct)/2.  Furthermore,  =  l/(C)l(^o/7’s)exp(-r2/r^),  where  the  spotsize  is 

given  by  =  ro(l  +  To  zeroth  order,  the  envelope  propagates  at  c,  i.e.,  group 

velocity  effects  are  not  included.  The  first  order  solution  is  =  /(C  -  i'tp')ap{r,  ko,T]), 
as  indicated  by  Eq.  (31),  where  hp  =  exp(V’)  and  V’  and  -ip'  are  given  by  Eq.  (26).  The 
first  order  solution  contains  group  velocity  effects,  with  the  group  velocity,  Eq.  (35),  in 
agreement  with  the  exact  expression,  Eq.  (19),  to  leading  order.  Higher  order  solutions  for 
a  contain  corrections  to  higher  order  in  the  small  parameter  l/k^L"^. 

Of  fundamental  importance  is  the  group  velocity  of  a  tightly  focused,  ultrashort  laser 
pulse.  Solutions  to  the  paraxial  wave  equation  imply  an  envelope  velocity  Ve{r,  z)  given  by 
Eq.  (7).  This  can  not  be  interpreted  as  a  group  velocity,  since  there  are  regions  in  {r,z) 
for  which  Ve  >  c.  At  the  focal  point,  1  -  Ve(0,0)/c  ~  2/fcgrg,  a  value  which  differs  firom 
1  -  Vg/c  by  a  factor  of  two.  By  analyzing  the  first  order  solution  d^^^  to  the  full  wave 
equation,  the  interpretation  of  the  paraxial  result  Ve(r,  z)  becomes  clear.  Using  d(^\  the 
expression  for  the  local  pulse  centroid  (C)c,  which  is  a  function  of  r  and  77  as  given  by 
Eq.  (33),  was  determined  by  averaging  C  over  the  axial  pulse  profile  at  a  fixed  radius  r. 
If  the  motion  of  the  local  centroid  (C)c,  which  is  a  function  of  r  and  77,  is  used  to  define 
a  local  envelope  velocity  7)e(r,  77),  the  result,  Eq.  (32),  is  identical  to  the  paraxial  result 
Ve{r,  z)  with  z  ^  77.  It  clear  that  Ue(r,  77)  is  not  the  pulse  group  velocity,  but  instead  is  the 
velocity  of  the  local  centroid  (C)c  at  a  fixed  radius.  The  fact  that  regions  in  (r,  77)  exist 
where  Ve{r,r))  >  c  indicates  that  Ue(r,  77)  does  not  represent  the  axial  velocity  at  which 
pulse  energy  is  transported. 

An  improved  definition  of  the  pulse  group  velocity  is  given  by  considering  the  motion 
of  the  centroid  C  of  the  entire  pulse,  Eqs.  (16)  and  (34),  which  is  determined  by  averaging 
C  over  both  the  axial  and  radial  pulse  profiles.  The  pulse  group  velocity  is  given  by 
Vg  =  dz/dt,  where  C  ^  z  -  ct.  The  exact  expression  for  Vg  is  given  by  Eq.  (19).  In  the 
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limits  fcorg  ^  1  and  :§>  1, 


^0^0  V  ^  ^0-^  / 


(42) 


In  addition  to  finite  spotsize  and  plasma  dispersive  effects,  the  group  velocity  is  also  slightly 
reduced  due  to  the  finite  pulse  length.  Notice  that  Vg  is  independent  of  2  such  that  the  pulse 
centroid  propagates  at  a  constant  velocity.  Thus  Eq.  (42)  is  valid  both  at  focus  and  far  from 
focus.  The  relativistic  factor  associated  with  the  pulse  group  velocity,  7s  =  (l-uJ/c^) 
to  leading  order,  is 

79  -  \/2(7rro/A)  (l  +  .  (43) 


For  a  tightly  focused  laser  pulse,  7g  can  be  relatively  small. 

As  an  example,  consider  the  possibihty  of  accelerating  an  electron  from  rest  with 
the  axial  ponderomotive  force  of  an  intense,  ultrashort  laser  pulse.  For  a  fixed  axial 
ponderomotive  force,  electron  trapping  is  easier  the  lower  the  value  of  jg.  For  a  laser  pulse 
in  vacuum  with  tq/A  =  3,  7^  ~  13  and  it  may  be  possible  for  a  sufliciently  intense  laser 
pulse  to  pick  up  and  accelerate  an  electron  from  rest. 
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Fig.  1:  Second  order  solution:  Real  part  of  a  =  dexpfifcoC)  plotted  versus  r  and  (  =  z—ct 
at  (a)  t  =  0  and  (b)  t  =  2Zjilc  for  L/X  =  tq/X  =  2  and  A  =  2.5  nm.  This  shows 
the  evolution  of  the  normahzed  vector  potential  of  the  laser  field  for  a  pulse  of  15 
fs  duration  with  a  central  wavelength  of  2.5  nm.  Plot  (a)  shows  the  initial  laser 
field  at  t  =  0  at  the  focal  point  where  the  minimum  spot  size  is  tq  =  5  /xm.  Plot 
(b)  shows  the  laser  field  after  propagating  a  distance  of  two  Rayleigh  lengths  from 
the  focal  point,  ct  =  2Zr  ~  63  /xm.  Notice  in  (b)  the  curvature  of  the  wavefronts 
as  the  pulse  diffracts. 
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Fig.  2:  Percentage  difference  100(a("')  -  between  numerical  solution 

to  the  wave  equation  and  n-th  order  analytic  solution  plotted  versus 

r  and  Q  -  z  -  ct  iox  (a)  the  first  order  solution  and  (b)  the  second  order 
solution  The  parameters  are  L/X  =  tq/X  =  2  and  A  =  2.5  uva.  and  the 

comparisons  are  done  after  propagating  a  distance  r;  =  (z  +  ct)/2  =  Zr  from  the 
focal  point. 


id  (b)  V  =  Ve,  where  Ve  =  Ve{r  =  0,rj)  is  the  local  envelope  velocity 
)]  evaluated  along  the  axis  r  =  0.  Here,  L/X  =  tq/X  =  2,  X  =  2.5 
3  the  focal  point. 


